Approximation Orders for Interpolation by Surface 
Splines to Rough Functions 



dJiI and Will Lighiy 



Rob Brownleco and Will Ligh 
Department of Mathematics and Computer Science, University of Leicester, University 

Road, Leicester LEI 7RH, England. 

Abstract 

In this paper we consider the approximation of functions by radial basic function inter- 
polants. There is a plethora of results about the asymptotic behaviour of the error between 
appropriately smooth functions and their interpolants, as the interpolation points fill out a 
bounded domain in M'^. In all of these cases, the analysis takes place in a natural function 
space dictated by the choice of radial basic function - the native space. In many cases, the 
native space contains functions possessing a certain amount of smoothness. We address the 
question of what can be said about these error estimates when the function being inter- 
polated fails to have the required smoothness. These are the rough functions of the title. 
We limit our discussion to surface splines, as an exemplar of a wider class of radial basic 
functions, because we feel our techniques are most easily seen and understood in this setting. 
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1 INTRODUCTION 



The process of interpolation by translates of a basic function is a popular tool for the 
reconstruction of a multivariate function from a scattered data set. The setup of the problem 
is as follows. We are supplied with a finite set of interpolation points A C and a function 
f : A^R. We wish to construct an interpolant to / of the form 

{Sf){x) = Y,f^ai^i^-a)+pix), foTxeR''. (1.1) 

Here, t/i is a real-valued function defined on R^, and the principle ingredient of our inter- 
polant is the use of the translates of ip by the points in A. The function tp is referred to as 
the basic function. The function p in Equation (|l.ip is a polynomial on H*^ of total degree 
at most k — 1. The linear space of all such polynomials will be denoted by Hk-i- Of course, 
for Sf to interpolate / the real numbers and the polynomial p must be chosen to satisfy 
the system 

{Sf){a) = /(a), for aeA. 

It is natural to desire a unique solution to the above system. However, with the present 
setup, there are less conditions available to determine Sf than there are free parameters 
in Sf. There is a standard way of determining the remaining conditions, which are often 
called the natural boundary conditions: 

^ /iag(a) = 0, for aU q G Hfe_i. 

aeA 

It is now essential that A is H/j_i-unisolvent. This means that if q £ Hk-i vanishes on A 
then q must be zero. Otherwise the polynomial term can be adjusted by any polynomial 
which is zero on A. However, more conditions are needed to ensure uniqueness of the 
interpolant. The requirement that -0 should be strictly conditionally positive definite of 
order k is one possible assumption. To see explanations of why these conditions arise, 
the reader is directed to Cheney & Light (2000). In most of the common applications the 
function -0 is a radial function. That is, there is a function (f> : IR4- — > IR such that 0; = 0o | ■ |, 
where | • | is the Euclidean norm. In these cases we refer to as a radial basic function. 
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Duchon (1976, 1978) was amongst the first to study interpolation problems of this flavour. 
His approach was to formulate the interpolation problem as a variational one. To do this 
we assume we have a space of continuous functions X which carries a seminorm | • | . The 
so-called minimal norm interpolant to f Cz X on A from X is the function Sf G X satisfying 

1. (5/) (a) = /(a), for aU a e A; 

2- \Sf \ < \g\, for all g e X such that g{a) = f{a) for aU a e A. 

The spaces that Duchon considers are in fact spaces of tempered distributions which he is 
able to embed in C(]R''). Let §' be the space of all tempered distributions on M''. The 
particular spaces of distributions that we will be concerned with are called Beppo-Levi 
spaces. The k*"^ order Beppo-Levi space is denoted by BL*'{Vl) and defined as 

BL^{n) = {/ e S' : D"/ e L2{n), a e TL^, \a\ = , 

with seminorm 

\f\k,n - f E / dx] , / G BL'in). 

The constants Ca are chosen so that the seminorm is rotationally invariant: 

J2 CaX^" = for all X G M'*. 

\a\—k 

We assume throughout the paper that 2k > d, because this has the affect that BL^{Vl) is 
embedded in the continuous functions (Duchon 1976). The spaces give rise to 

minimal norm interpolants which are exactly of the form given in Equation (|l.ip . where the 
radial basic function is a; jxp*''^'' or a; i-^ log |a;|, depending on the parity of d. 

It is perhaps no surprise to learn that the related functions tp are strictly conditionally 
positive definite of some appropriate order. The name given to interpolants employing these 
basic functions is surface splines. This is because they are a genuine multivariate analogue 
of the well-loved natural splines in one dimension. 

It is of central importance to understand the behaviour of the error between a function 
/ : ^ H and its interpolant as the set A (1 ^ becomes "dense" in Q,. The measure of 
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density we employ is the fill-distance h = sup^.^^ minag^|a; — a\. One might hope that for 
some suitable norm || ■ || there is a constant 7, independent of / and h, such that 

\\f-Sf\\ = 0{h^), as/i^O. 

In the case of the Beppo-Levi spaces, there is a considerable freedom of choice for the norm 
in which the error between / and Sf is measured. The most widely quoted result concerns 
the norm || • || 1,^(0), but for variety we prefer to deal with the Lp-norm. To do this it is 
helpful to assume is a bounded domain, whose boundary is sufficiently smooth. In this 
case there is a constant C > 0, independent of / and h, such that for all / G BL''{Q), 

There has been considerable interest recently in the following very natural question. 
What happens if the function / does not possess sufficient smoothness to lie in BL''{fl)l It 
may well be that / lies in BL"^{^1), where 2k > 2m > d. The condition 2m > d ensures 
that /(a) exists for each a £ A, and so Sf certainly exists. However, \f\k,n is not defined. 
It is simple to conjecture that the new error estimate should be 

ii/-s/iiM.,<{s:|/t*Jf'"* ?i^;r- ('^^) 

It is perhaps surprising to the uninitiated reader that this estimate is not true even with 
the reasonable restrictions we have placed on k and m. We are going to describe a recent 
result from Johnson (2002). To do that, we recall the familiar definition of a Sobolev space. 
Let W2^(n) denote the fc*^ order Sobolev space, which consists of functions all of whose 
derivatives up to and including order k are in L2{Cl). It is a Banach space under the norm 

ll/IU,a = ( El/I'n 1 , where / e W^{n). 

We have already tacitly alluded to the Sobolev embedding theorem which states that when 
O is reasonably regular (for example, when Cl possesses a Lipschitz continuous boundary) 
and k > d/2, then the space (O) can be embedded in C(0) (see Adams 1978, Theorem 
5.4, p. 97). Now Johnson's result is as follows. 



Theorem 1.1 (Johnson). Let he the unit ball in H*^ and assume d/2 < m < k. For every 
ho > 0, there exists an / G W^iJR'^) and a sequence of sets {Anjnen with the following 
properties: 

(i) each set An consists of finitely many points contained in fJ; 

(ii) the fill-distance of each set An is at most ho; 

(Hi) if S'^f is the surface spline interpolant to f from BL^iJM^) associated with An, 
for each n e IN, then \\S]tf\\Li{n) — > oo as n ^ oo. 

If the surface sphne interpolation operator is unbounded, there is of course no possibihty 
of getting an error estimate of the kind we conjectured. Johnson's proof uses point sets which 
have a special feature. We define the separation distance of An as qn = min{|a — h\/2 : a,b G 
An, a ^ b}. Let the fill-distance of each An be /i„. In Johnson's proof, the construction of An 
is such that qn/hn — > 0. We make this remark, because Johnson's result in one dimension 
refers to interpolation by natural splines, and in this setting the connection between the 
separation distance and the unboundedness of SJ! has been known for some time. What 
is also known in the one-dimensional case is that if the separation distance is tied to the 
fill-distance, then a result of the type we are seeking is true. Theorem 13.51 is the definitive 
result we obtain, and is the formalisation of the conjectured bounds in Equation (|1.3p . 

Subsequent to carrying out this work, we became aware of independent work by Yoon 
(2002). In that paper, error bounds for the case we consider here are also offered. Because 
of Yoon's technique of proof, which is considerably different to our own, he obtains error 
bounds for functions / with the additional restriction that / lies in W^{il), so the results 
here have wider applicability. However, Yoon does consider the shifted surface splines, whilst 
in this paper we have chosen to consider only surface splines as an exemplar of what can 
be achieved. At the end of Section 3 we offer some comments on the difference between our 
approach and that of Yoon. 

To close this section we introduce some notation that will be employed throughout the 
paper. The support of a function (f> : fR"* ^ IR is defined to be the closure of the set 
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{x G II : (j){x) ^ 0}, and is denoted by supp((/)). The volume of a bounded set is the 
quantity J^^ dx and wiU be denoted vol(ri). We make much use of the space Ilm-i, so for 
brevity we fix I as the dimension of this space. FinaUy, when we write / we mean the Fourier 
transform of /. The context wiU clarify whether the Fourier transform is the natural one 
on Li(]R'^): 

or one of its several extensions to i2(Il'') or §'. 

2 SOBOLEV EXTENSION THEORY 

In this section we intend to collect together a number of useful results, chiefly about the sorts 
of extensions which can be carried out on Sobolev spaces. We begin with the well-known 
result which can be found in many of the standard texts. Of course, the precise nature of 
the set i7 in the following theorem varies from book to book, and we have not striven here 
for the utmost generality, because that is not really a part of our agenda in this paper. 

Theorem 2.1 (Adams 1978, Theorem 4.32, p. 91). Let he an open, bounded subset of 
H'* satisfying the uniform cone condition. For every f G W2^{fl) there is an G ^^"(IR'*) 
satisfying \n= f . Moreover, there is a positive constant K = K[il) such that for all 

f e w^ifi), 

Wf^Wni.K' ^ K\\f\\m,n- 

We remark that the extension can be chosen to be supported on any compact subset 
of M'* containing Q. To see this, we construct in accordance with Theorem 12.11 then 
select rj G C™(]R'') such that rj{x) ~ 1 for x G fl. Now, if we consider the compactly 
supported function /q^ = r;/^ G ^^"(IR''), we have /o'lij^ /. An elementary application of 
the Leibniz formula gives 

ll/o''lL,R'' < C||/||™,n, where C = C{n,ij). 

One of the nice features of the above extension is that the behaviour of the constant 
K{n) can be understood for simple choices of il. The reason for this is of course the choice 
of Q and the way the seminorms defining the Sobolev norms behave under dilations and 
translations of f2. 
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Lemma 2.2. Let ^ be a measurable subset o/H'*. Define the mapping a : R** — > H"^ by 
a{x) = a + h{x -t), where h>0, and a, t, x G'R'''. Then for all f e W^{a{Q)), 

Proof. We have, for |a| = m, 

{D'-ifoa)){x) = h"'iD'^f){a{x)). 

Thus, 

l/°<n= E / \iD"{foa)){x)\'dx 

\a\=m 

= h'^ c^l \{D'^f){cj{x))\^dx. 

I I JQ. 
\a\—m 

Now, using the change of variables y = a{x), 

l/°^l™.o = /^" E / \iD"f)iy)\'dy = h'-^-VtMny □ 
|a|=™ -^-C") 

Unfortunately, the Sobolev extension refers to the Sobolev norm. We want to work with 
a norm which is more convenient for our purposes. This norm is in fact equivalent to the 
Sobolev norm, as we shall now see. 

Lemma 2.3. Let Ct be an open subset of R"^ having the cone property and a Lipschitz- 
continuous boundary. Let bi, . . . ,be G Cl be unisolvent with respect to Hm-i- Define a norm 
on W^{n) via 

/ I \ 1/2 

ll/lln = f + El/(^^)l' ) ' / e W^(Sl). 

There are positive constants Ki and K2 such that for all f e W2^{Q), 

i^i||/IU,n< ||/||o<i^2||/|U,n. 

Proof. The conditions imposed on m and fl ensure that W2^{fl) is continuously embedded 
in C{Cl) (Adams 1978, Theorem 5.4, p. 97). So, given x £ Cl, there is a constant C such 
that |/(a;)| < C'\\f\\m,n for all / £ W^{fl). Thus, there are constants Ci, . . . ,Ce such that 

11/11^ < \ft,a + E Ci\\ffm,n < (1 + E ^0 11^"™.^- ^2-^) 

1=1 i=l 
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On the other hand, suppose there is no positive number K with ||/||m,si < -^||/||si for aU 
/ e W^{^). Then there is a sequence {/j} in W^'{n) with 

||/,||™,a = l and ||/,||a < ^, for j = 1, 2, . . . . 

The ReUich selection theorem (Braess 1997, Theorem 1.9, p. 32) states that W™{Q) is 
compactly embedded in W^~^ [Vl) . Therefore, as {/j} is bounded in W^{^), this sequence 
must contain a convergent subsequence in W^~^{Vl). With no loss of generality we shall 
assume {/j} itself converges in W^~^{p). Thus {/j} is a Cauchy sequence in W^™~^(ri). 
Next, as ||/j||o — * it follows that |/j|m.o — » 0. Moreover, 

ll/j ^ /fcllm,n — ll/j ^ /fc|lm-l,fl + l/j ^ /fe|m,0 

Since {/,} is a Cauchy sequence in W^™~^(il), and |/j|m,si ^ 0, it follows that {/j} is a 
Cauchy sequence in W2^{VL). Since VK™(ri) is complete with respect to || • ||m,n, this sequence 
converges to a hmit / e W™(f2). By Equation p.l|l . 

and hence ||/ — fj\\n as j — * oo. Since ||/j||o ^ 0, it follows that / = 0. Because 
||/il|m,n = 1- J = li2,..., it follows that ||/||m.o = 1- This contradiction establishes the 
result. □ 
We are almost ready to state the key result which we will employ in our later proofs about 
error estimates. Before we do this, let us make a simple observation. Look at the unisolvent 
points bi, . . . ,b( in the statement of the previous Lemma. Since W^{Q) can be embedded 
in C{fl), it makes sense to talk about the interpolation projection P : W^{il) — > Ilm-i 
based on these points. Furthermore, under certain nice conditions (for example f2 being a 
bounded domain), P is the orthogonal projection of W^{il) onto 11^-1. 

Lemma 2.4. Let B be any ball of radius h and center a £ H'', and let f £ W:^{B). 
Whenever bi,. . . ,be G are unisolvent with respect to Hm-i let Pb ■ C(IR'') n„,_i be 
the Lagrange interpolation operator on 6i, . . . , 6^. Then there exists c — (ci, . . . , q) G 
and g e W^(R'^) such that 
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1. g{x) = {f -Pcf){x) for allxeB; 

2. g{x) = for all |x - a| > 2h; 

3. there exists aC > Q, independent of f and B, such that \g\m,'R'^ ^ C|/|m,s- 
Furthermore, Ci, . . . , q can he arranged so that Ci = a. 

Proof. Let Bi be the unit ball in M!^ and let B2 = 2Bi . Let hi, . . . ,bi £ i?i be unisolvent 
with respect to lim-i- Define cf{x) = h~^{x — a) for all x £ H''. Set q = a^^{hi) for 
i = 1, . . . ,i so that ci, . . . ,ce & B are unisolvent with respect to n.,„-i. Take / G W2''{B). 
Then (/ - Pcf) o CT-i e VF2"(Bi). Set F = {f - Pcf) o cr-i. Let F^^ be constructed as an 
extension to F on Bi. By Theorem 12. II and the remarks following it, we can assume F^^ is 
supported on B2. Define g ^ F^^ o a e W:^'{R'^). Let a; e B. Since (j{B) = Bi there is a 
y & Bi such that x — a^^{y). Then, 

g{x) = [F^- o a) (a;) = F^-{y) = ((/ - PJ) o a-^){y) = (/ - 

Also, for X e H'^ with jx — a| > 2h, we have |cr(x)| > 2. Since is supported on B2, 
g{x) = for |x — a| > 2h. Hence, g satisfies properties 1 and 2. By Theorem 12. II there is a 
Ki, independent of / and B, such that 

We have seen in Lemma [2731 that if we endow W^{Bi) and W^{B2) with the norms 

/ e \ 1/2 

then [| • and || • \\rn.Bi sue equivalent for i = 1,2. Thus, there are constants K2 and K3, 
independent of / and B, such that 

< K2\\F''^\\m,B, < KiK2\\Fln,B, < i^lifsi^S Ili^IUi • 

Set C = K1K2K3. Since F^^{h,) = F{h,) = (/ - PJ){a-\b,)) = (/ - = for z = 

it follows that \F^'\rn.B, < C\F\,n,B,- Thus, \goa-^L,M^ < C\{f^PJ)oa-^\„r,B,- 
Now, Lemma |2. 21 can be employed twice to give 

\gL,u-=h''/'-"'\goa-'\„^^^,<Ch^/'-'-\{f-PJ)oa-'UB,=C^ 
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Finally, we observe that \ f — Pcf\m,B — \ f\m,B to complete the first part of the proof. The 
remaining part follows by selecting bi — Q and choosing b2, ■ ■ ■ ,bi accordingly in the above 
construction. □ 

Lemma 2.5 (Duchon 1978). Let ft be an open, bounded, connected subset of having 
the cone property and a Lip schitz- continuous boundary. Let f G W^{fl). Then there 
exists a unique element f^^ £ BL"^(JR'^) such that f , and amongst all elements 

of BL'^{K^) satisfying this condition, \f^\mM'^ minimal. Furthermore, there exists a 
constant K = K{Vl) such that, for all f G W^{^), 

3 ERROR ESTIMATES 

We arrive now at our main section, in which we derive the required error estimates. Our 
strategy is simple. We begin with a function / in BL™{Ji'^). We want to estimate ||/ — 5*^/11 
for some suitable norm || • ||, where Sk is the minimal norm interpolation operator from 
BL^{K^), and fc > m. We suppose that we already have an error bound using the norm 
II • II for all functions g £ BL''(R'^). Our proof now proceeds as follows. Firstly, we adjust 
/ in a somewhat delicate manner, obtaining a function F, still in i?L™(]R'^), and with 
semmorm m BL™(]R'^) not too far away from that of /. We then smooth F by convolving 
it with a function e C^(]R''). The key feature of the adjustment of / to F is that 
(0 * F){a) — f{a) for every point a in our set of interpolation points. It then follows that 
F e BL'^iM.'^). We then use the usual error estimate in BL''(R'^). A standard procedure 
(Lemma [331 then takes us back to an error estimate in i3L''"(Il'^). 

Lemma 3.1. Let m < k and let (p G C^{Ji'^). For each h > 0, let (j)h{x) = h~'^(j){x/h) 
for X G . Then there exists a constant C > 0, independent of h, such that for all 
f G BL"'(R'^), 

Furthermore, we have \(j)h * f\k tr'' — o(/i™^'^) as h 0. 
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Proof. The chain rule for differentiation gives {D^(l)h){x) = h''^'^+^^\\D^(l)){x/h) for all 
X e R'', and 7 £ Thus, for [3 e7L'\_ with |/31 = m we have 



\{D'^<i}h*DPf){x)\^ dx^ 



/ (D^0,O(^-y)P''/)(y) 



/j-2(d+|7l) 



/,-2|7l 



^/,-2|7l 



{D-^ (t)){t){DP f){x - ht) dt 

d 

[D'' (j)){t)[D'^ f){x - ht) dt 



dx 



dx 



dx, (3.1) 



where K = supp (0). An application of the Cauchy-Schwartz inequality gives 



\{D''^h * D^f){x)\^ dx < /i-2|7l / / |(Z?70)(i)|2 ^Al \^D^ f)^x - ht)\^ dt dx, 

E," JW' \ Jk J \ Jk I 



and so, 



\{D-^(l)h*D'^f){x)\'^ dx <h-'^^-^^ / / / \{D'^f){x-ht)\'^dtdx. (3.2) 



The Parseval formula together with the relation {D°'{<j)h * f ))^ = {i ■ )°^{4>h * /)^ provide us 
with the equality 



V c„ / |(D«(</.„ * /))(x)p dx=y^c^ [ \{ixr[<PH * /)^(a:)P dx 
^R"^ ^R" 

= / c^x^''Uh*fr{x)\^ dx. (3.3) 



\a\=k 



\a\ = k 



Now, when Equation p.3p is used in conjunction with the relation 



E2a I |2fc I |2(m+fe-m) 2/3 
CaX =\x\ = |a;| ^ I = cpx ' c^x 
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\a\=k 



m I7I— fc— m 
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we obtain 



I I 7 '^H'^ -^H^ 1^1 I I f 

^ / ^ c,x2'^i(zx)''(</>,*/r(x)pdx 



\p\—m \^\—k—m 



= E E cJ \{ixr{Dp{4>H*f)r{x)\' dx 

= E E cJ \{D'<{D'\cpn*f))r{^)?dx 

\[3\=m |7|=fe-m 

= E E cJ \{D'<{Df\cpn*fmx)?dx. 

|/3|=m |7|=fe-m "^"^'^ 

Since the operation of differentiation commutes with convolution, we have that 

CcJ J{D''{4>h*mx)\^ dx^ E J J{D^4>h*D^f){x)\^ dx. (3.4) 

\a\ — k ^ m \-y\—k—m ^ 

Combining Equation (|3.2p with Equation (|3.4p we deduce that 

< V V C7/i^2|7l /■ \{^D^cp){t)\^ dt I [ \{D^ f){x - ht)\^ dtdx 



\a\—k 



\f3\—rn \^\—k—rn 



Fubini's theorem permits us to change the order of integration in the previous inequahty. 
Thus, 

Y^o^f \{D%cl^H*mxtdx<e^^-^Ml-^^^Y^pl I \{D^f){x~ht)\'dxdt. 
I I I, J^"^ ' ifli JkJr'' 

\a\—k \i3\—m 

FinaUy, a change of variables in the inner integral above yields 

E^a/ \{D'^{cl^H^mx)\'dx<h^("^-'^'^\cj^\l_^^, Y [ I \{D^m)?dzdt. 

\a\—k \f3\—m 

Setting C = |(/)|j._,„ j^d •\/vo1(k7 we conclude that \4ih*f\k,w' — '^l/lm,]R'' as required. 

To deal with the remaining statement of the lemma, we observe that for 7 7^ we have 

(D-'^m dt^ f (D-'^m dt = {wcj^m = {{z ■ r^m = o. 
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Then it follows from Equation (|3.ip that for |/3| — m, 

\{D''^h*D^f){x)\^ dx = h-'^'>\ f f {D'>^){t){{D^f){x-ht)-{D^f){x)) dt dx. 
Now, if we continue in precisely the same manner as before, we obtain 

<;^2(™-.)|^|2_^^^ ^ ^ / / \iD''f)ix-ht)-{D^f){x)\' dxdt. 

Since L"^/ G L2(]R^) for each /? e S| with |/3| = m, it follows that for almost aU t,x e IR'', 
\{D'^f){x~ht)^{D'^f){x)\^0, as/i">0. 

Furthermore, setting 

g{x,t) = 2\{D^f){x - + 2|(D^/)(a;)p, for almost aU x,t G M'*, 
we see that 

|P'3^)(x-/it)-(i?^/)(x)|2<5(x,t), 
for almost all x.t G H'^ and each h > 0. It follows by calculations similar to those used 



above that 

/ / g{x,t)dxdt^4vo\{K) \{D'^ f){x)\^ dx < oo. 
Applying Lebesgue's dominated convergence theorem, we obtain 

\{D'^f){x - ht) - {Dl^f){x)\'^ dxdt ^0, as /i ^ 0. 



IK JR<* 

Hence, for m < k, \(f)h * f\k — o(/i™~'^) as /i — » 0. □ 
Lemma 3.2. Suppose (f> e C^(R^) is supported on the unit ball and satisfies 



4>{x) dx — 1 and I (j){x)x°' dx — 0, for all < |a| < m — 1. 

For each e > and x G R'^, let (j)eix) = e-'^(l){x/e). Let B he any ball of radius h and center 
a G ■ For a fixed p G Ilm-i let f be a mapping from IR'^ to IR such that f{x) — p{x) for 
all X £ B. Then {(pe * f){a) — p{a) for all s < h. 
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Proof. Let Bi denote the unit ball in M''. We begin by employing a change of variables to 
deduce 

(0£ * /)(a) = / (l)e{a~y)f{v) dy 

4>{^)f{y)dy 

'x)f{a — xe) dx 
x)f{a — xe) dx. 

Then, for x G |(a — xe) — a| < e < h. Thus, /(a — xe) = p{a — xe) for all x G Bi. 
Moreover, there are numbers 6q such that p{a — xe) — p{a) + X]o<|a|<m-i ^ax" . Hence, 

{4>e * f){a) = / (j){x)p{a - xe) dx 



= f (j){x) (p{a) + V b„xA dx 

= pia). □ 
Definition 3.3. Let be an open, bounded subset o/R'^. Let A be a set of points in Q. 
The quantity sup^^jj inf(jg_4|a; — a\ = h is called the fill-distance of A in f2. The separation 
of A is given by the quantity 



q — mm 

a.beA 2 

The quantity h/q will be called the mesh-ratio of A. 

Theorem 3.4. Let A be a finite subset ofM!^ of separation q > and let d < 2m < 2k. 
Then for all f £ BL™(]R'^) there exists an F e BL''(R'^) such that 

L F{a) = f{a) for all a e A; 

2. there exists a C > 0, independent of f and q, such that < C'|/|„ j^d and 

Proof Take / £ SL'"(IR''). For each ae A\et BaC denote the ball of radius S = q/A 
centered at a. For each Ba let ga be constructed in accordance with Lemma [2.41 That is, 
for each a e A take c' = (c2, . . . , q) G B^~^ and ga G W^2™(I^'') such that 

1. a, C2, . . . , ce are unisolvent with respect to n„j_i; 
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2. ga{x) = (/ - P(a,c')f){x) for all X £ Ba\ 

3. P(a.c')f e n,„_i and {P(a,c')f){a) = /(a); 

4. = for all |a: - a| > 26; 

5. there exists a Ci > 0, independent of / and Ba, such that \ga\m e/* — Cil/lm.Sa- 
Note that if a ^ b, then supp (ga) does not intersect supp (gt), because if a; G supp (ga) then 

\x-b\> \b- a\-\x~ a\>2q-25 ^ 66. 
Using the observation above regarding the supports of the g^s it follows that 



dx 



supp (gb) 



J2iD"ga){x) 



a£A 



dx 



supp (Sb) 



a|=m b£A 

= E ^"E 

— Igalm.R'^- 

Applying Condition 5 to the above equality we have 

2 

ESa <CfY.\f\l,s^<Cf\f\l,,,. 
aeA "i.R'' ae^l 

Now set = / — It then follows from Condition 2 above that i/(a;) — {P(a.c')f){x) 
for all X £ Ba, and from Condition 3 that H{a) = /(a) for aU a e A. Let cj) e C(f (M'') be 
supported on the unit ball and enjoy the properties 



/ dx = 1 and / 4){x)x°' dx 

Jtr.'' Jw 



= 0, for all < lal < m - 1. 



Now set F — (j>s * H . Using Lemma [3. 1[ there is a constant C2 > 0, independent of q and 
/, such that 



/ - E^a 



aeA 



^ a£A 

<2C2{l + C!)6'("'-'^\f\l^,. 
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Similarly, there is a constant C3 > 0, independent of q and /, such that 

2 



f -^9a 



<2C^{1 + Cl)\f\l^^,. 

Thus li^lfc^Rd < C'?'"~''|/lm,E,'' ^iid l-F'Im.E,'' < C\f\m,M'^ some appropriate constant 
C > 0. Since F = (ps * H and -ff |_b^ G Hm-i for each a G it follows from Lemma [5T^ that 
F{a) = H{a) = /(a) for all a e A. D 

Theorem 3.5. Let D, be an open, bounded, connected subset of H'^ satisfying the cone 
property and having a Lipschitz-continuous boundary. Suppose also d < 2m < 2k. For each 
h > 0, let Ah be a finite, Ilk-i-unisolvent subset of Q with fill-distance h. Assume also 
that there is a quantity p > such that the mesh-ratio of each Ah is bounded by p for all 
h > 0. For each mapping f : Ah 11, let S^f be the minimal norm interpolant to f on 
Ah from BL^iR^). Then there exists a constant C > 0, independent of h, such that for all 

f e BL^{n), 

Proof. Take / G BL'^in). By Duchon (1976), / G W:^"'{n). We define in accordance 
with Lemma 12.51 For most of this proof we wish to work with and not /, so for 
convenience we shall write / instead of Construct F in accordance with Theorem 13.41 
and set G = f — F. Then F{a) = /(a) and G{a) — for all a G Ah- Furthermore, there is 
a constant Ci > 0, independent of / and h, such that 

/ h \ 

\F\k,n'' < C^l {-] l/lm.R'': (3-5) 

I^L^Kd < |/|„i,]R,<i + \P\m,TR'' < (1 + (3.6) 

Thus S'l^f — S'l^F and S^'^^G — 0, where we have adopted the obvious notation for S*^. 
Hence, 

ll/-^.VllL,(a) = ll/-^fc'i^llL,(o) - \\F+G-S'^F\\L^^n) < \\F-S'^F\\L^^a) + \\G-StGh^^n)- 
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Now, employing Duchon's (1978) error estimates for surface splines p.2p . there are positive 
constants C2 > and C3 > 0, independent of h and /, such that 

where we have defined 



Pij) = 



J-| + |, 2<p<oo 



J, l<p<2 ■ 

Finally, using the bounds in Equations p. 51) and p. 61) we have 

11/ - 5[7IIl,(o) < C4/i'='(™)|/U,R., ash^O, 

for some appropriate C4 > 0. To complete the proof we remind ourselves that we have 
substituted with /, and so an application of Lemma [2.51 shows that we can find C5 > 
such that 

11/ - ^[71Il,(o) < C,hP^"'^n^^^. < C,C,h^^"^^\fU,n, as /i ^ 0. □ 

We conclude this section with a brief commentary on the approach of Yoon (2002). It 
is hardly surprising that Yoon's technique also utilises a smoothing via convolution with a 
smooth kernel function corresponding closely to our function used in the proof of The- 
orem 13.51 However, Yoon's approach is simply to smooth at this stage, obtaining the 
equivalent of our function F in the proof of Theorem 13. 51 Because there is no preprocessing 
of / to H, Yoon's function F does not enjoy the nice property F{a) = /(a) for all a E A. It 
is this property which makes the following step, where we treat G ^ f — F, a fairly simple 
process. Correspondingly, Yoon has considerably more difficulty treating his function G. 
Our method also yields the same bound as that in Yoon, but for a wider class of functions. 
Indeed we would suggest that BL"^{Vl) is the natural class of functions for which one would 
wish an error estimate of the type given in Theorem 13.51 
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